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The effect of an external bias voltage and fluctuating electromag- 
netic fields on both the fundamental frequency and damping of can- 
tilever vibrations is considered. An external voltage induces surface 
charges causing cantilever-sample electrostatic attraction. A similar ef- 
fect arises from charged defects in dielectrics that cause spatial fluctu- 
ations of electrostatic fields. The cantilever motion results in charge 
displacements giving rise to Joule losses and damping. It is shown that 
the dissipation increases with decreasing conductivity and thickness of 
the substrate, a result that is potentially useful for sample diagnostics. 
Fluctuating electromagnetic fields between the two surfaces also induce 
attractive (Casimir) forces. It is shown that the shift in the cantilever 
fundamental frequency due to the Casimir force is close to the shift ob- 
served in recent experiments of Stipe et al. [1]. Both the electrostatic 
and Casimir forces have a strong effect on the cantilever eigenfrequen- 
cies, and both effects depend on the geometry of the cantilever tip. We 
consider cylindrical, spherical, and ellipsoidal tips moving parallel to 
a flat sample surface. The dependence of the cantilever effective mass 
and vibrational frequencies on the geometry of the tip is studied both 
numerically and analytically. 



I. INTRODUCTION 



The description of the motion of a microscopically small body near a massive 
solid is an important problem in various areas of modern physics, including atomic, 
scanning tunneling, and magnetic resonance force microscopies. The interactions be- 
tween bodies separated by a vacuum arise from electromagnetic fields of various types, 
the simplest example being an external potential difference applied between metallic 
bodies. In this case the surfaces in close proximity acquire opposite electric charges 
and experience an electrostatic attraction. A static electric field between different 
surfaces may exist even without any externally applied voltage; a zero-bias electro- 
static attraction may be caused by local variations in the work function or surface 
contamination resulting in an inhomogeneous electric field. This is referred to as the 
patch effect (see, for instance, References [2]- [4]). Also, spatial fluctuations of the 
electric field may be caused by charged defects in the bulk of a cantilever or a sample. 
In addition to these electrostatic effects, electrically neutral bodies interact via fluc- 
tuating electromagnetic fields. Fluctuating electric currents and charge densities (or 
polarizations) of a given body produce fields that act on other bodies. The strength 
of the interaction depends on the power spectrum of the fluctuations, which in turn 
are determined by the complex dielectric permittivity. Lifshitz treated this van der 
Waals force in a classic paper [5]. In the case of high-conductivity metal surfaces 
the effect is referred to as the Casimir force, which has been the subject of several 
recent experiments [6]- [9] and various surveys (see, for instance, Reference [10] and 
monographs [11], [12]). 

When a vibrating cantilever approaches a sample there is an increase in the dissi- 
pative force it experiences. Because of its great practical importance, the problem of 
noncontact friction has become a topic of increasing interest (see, for instance, Refer- 
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ences [13]- [20]). The utilization of micromechanical cantilevers could conceivably be 
extended considerably if the physical mechanisms responsible for the damping were 
understood. Such an understanding might lead to applications in which both the 
conservative (elastic or electrostatic) forces and the dissipative (velocity-dependent) 
forces are probed in order to obtain useful information about a sample. 

Attempts to explain the observed friction in terms of the van der Waals interaction 
[18], [21]- [23] have not met with much success [15], [19], [20]. On the other hand the 
effect of electric fields of various origins on the noncontact friction was shown in the 
experiments of Stipe et al. [1] to play an important role in determining the friction, 
although a microscopic theory still does not exist. It follows from the work of Stipe 
et al. [1] that any changes in the electric field (due, for example, to changes of tip- 
surface distance or variations in bias voltage associated with charge defects generated 
in dielectric substrates) produce significant modifications of the noncontact friction. 

In what follows we consider the perpendicularly oriented cantilever, which has 
been utilized as an ultrasensitive device for the study of forces and dissipation at 
small length scales [1]. The perpendicular orientation allows a much closer approach 
of the tip to the sample without jump-to-contact, and also much less cantilever damp- 
ing than in the parallel configuration. The influence of tip-substrate attraction on 
cantilever eigenfrequencies and Joule dissipation will be studied, and the dependence 
of these quantities on the tip geometry will be described both analytically and nu- 
merically. 

In Section II we develop expressions for the electric field , induced currents, and 
Joule heating in the case of a cylindrical tip. We consider both infinitely thick and 
finite-thickness substrates, and find that Joule heating is increased in the case of a 
thin substrate. Section III extends this analysis to the more realistic case of spherical 
or ellipsoidal tips, and in Section IV we consider the Casimir force, which can be the 
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dominant force at small tip-sample separations, for both cylindrical and ellipsoidal 
tips. In Section V we consider the force arising from spatial fluctuations of electro- 
static charge, and compare this force to the Casimir force for typical parameters of 
interest. Section VI is concerned with the effects of these attractive forces on the 
vibrational frequency and effective mass of the cantilever, and numerical results for 
these effects are given for typical values of various parameters. Our conclusions are 
summarized in Section VII. 



II. ELECTRIC FIELD AND INDUCED CURRENTS FOR A 
CYLINDRICAL TIP GEOMETRY 

In this section we consider first a model in which the cantilever tip is assumed 
to be a metallic cylinder with its axis parallel to a flat sample surface. The tip 
displacements are assumed to be nearly parallel to the surface, which will be the 
case when the oscillation amplitudes are sufficiently small. Figure 1 shows the model 
under consideration. The cantilever width b is taken to be much larger than the 
thickness a {b ^> a). Figure 1 shows a section in the plane of the cantilever motion. 
It is straightforward to obtain the static electric field distribution in the practically 
important case of small distances d such that the electrostatic field of the entire 
cylinder is effectively the same as that due only to its bottom part. (The criterion 
that d must satisfy for this to be the case is given below.) The problem is then 
reduced to solving the 2D Laplace equation with the boundary conditions that the 
potential has constant values 4>i an d 02 at the metallic surfaces. The electric field 
distribution outside the conductors is equal to the field due to a charge q and its 
image — q placed at a distance d\ = [(R + d) 2 — R 2 ] 1 / 2 from the surface. (Details may 
be found in Reference [24].) q represents charge per unit length along the y direction, 
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and is determined by the potential difference 4>i ~ 02 and the capacitance C given by 
C- 1 = 2cosh" 1 (f +d/R) = 2\n[(d + R + d 1 )/R] . (f) 

Thus [24] 

q = C(0i - 2 ) . (2) 



The electric field at a point r exterior to the tip and sample is given by 

- ri r - r 2 



E(r) = 2q 



|r — ri\ 2 |r — r 2 | 2 



(3) 



where ri and r 2 denote the positions of +q and -q, respectively. For the coordinate 
system shown in Figure 1, r 12 = ±d±e z , where e 2 is the unit vector along the z axis. 
The attractive cantilever-surface force can be calculated straightforwardly using Eqs. 
(l)-(3). We obtain 

T = bq 2 jd x (4) 

for the force on the tip. T is proportional to the length b of the cylindrical tip in our 
model, and is given by 

bV 2 R l l 2 

T = 2 7/2 # /2 > V = & - ^ ( 5 ) 

in the practically important case d <C R. 

Metals are almost inpenetrable by the electric field, the penetration being re- 
stricted to a very thin surface layer (Debye layer) of thickness R D considerably less 
than other characteristic lengths of the system, i.e., Rd <^ R,d,d\. Integrating the 
Poisson equation over a thin layer of the sample (but with thickness greater than 
Rd), we obtain the surface charge density 



where we have taken into account the boundary conditions E z (z = + ) = —4qdi/(x 2 + 
d\) and E z (z = 0") = 0. The last condition, of course, means that the overall field 
due to the charged tip and the surface (Debye) layer is zero in the volume of the 
metal substrate. 




FIG. 1. Cantilever-sample geometry for the model in which the cantilever "tip" is taken 
to be part of a cylinder of radius R. 

A somewhat different picture applies in the case of a moving charged tip. To 
be specific, we consider the case of oscillatory motion of the tip: X(t) = X cos uot. 
(Note that an oscillating cantilever is used, for example, for spin detection [25].) The 
cantilever charge is not changed when its tip moves parallel to the surface, while the 
sample charge varies in time at any fixed point. The varying sample charge induces 
an electric current, and the charge and current are related by the continuity equation, 
dp/dt + V • J = 0, where p and J are the charge and current densities, respectively. 
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Integration of the continuity equation over a small element of surface layer gives the 
relation between the normal current near the Debye layer and the time derivative of 
the surface charge density: 

d 

—a(x,t) = J z (x,z = ,t) . (7) 

We will assume that the cantilever motion is slow in the sense that oo <C 4irfi, where 
ji is the sample conductivity; in this case an iterative procedure described by Boyer 
[26] is applicable. The assumption of slow cantilever motion means that the surface 
charge follows the tip displacements practically instantaneously. To lowest order in 
the parameter o;/(47r/x), the surface charge may be approximated by Eq. (6) with the 
coordinate x replaced by x — X{t). This approximation corresponds to an effectively 
infinite value of volume conductivity and instantaneous recharging. Then, from Eq. 
(7), 



J z (x,z = , t) — /iE z (x, z — , t) = — — qdiuX sino>t— — ■ 2 . (8) 



The additional surface charge a' (up to first order in oo/mu) caused by the normal 
current (8) can easily be obtained by considering this charge as a source of the field 
E z (x,z = 0~,t). Taking into account the fact that the field E z induced by the 
additional surface charge has opposite signs below and above the surface z = 0, we 
obtain the perturbation of the surface charge density as follows: 

// ,x _ qd^Xpsmut d 1 , , 

a ^ t >~ 2 7T 2 /i dx (x - X f + d\- W 

For small tip oscillation amplitudes (X <C d±) we can ignore X(t) in the denominator. 
Thus we obtain the electric field due to the surface charge a': 

/oo Q 
dx'a'(x',t)—\n\r - e x x'\ . (10) 
oo OY 

Explicitly, the expression for the electric field that is first-order in u/fx is 



E'(r, t) 



quiX sin out d 



e x (\z\ + di) ± e z x 
x 2 + (\z\+d 1 ) 2 



) 



(11) 



where the + and - signs correspond to z > and z < 0, respectively. It follows from 
Eq. (11) that the characteristic length of the electric field variation is d±. Hence our 
simplifying assumption that the field formed by the full cylinder can be represented 
by only its bottom part (as in Figure 1) is valid in the case d± <C R. 

Equation (11) can be used to obtain the cantilever damping due to Joule dissipa- 
tion. The Joule power Wj is given by the integral 



Note that there is no dependence of Wj on the tip radius R. This results from 
the effects of two competitive factors determining the Joule loss. The first is the 
magnitude of the electric field in the bulk of the substrate, which increases with a 
decrease in the characteristic inhomogeneity size d 1 « (2dR) 1 ' 2 as R/d\. The second 
is the characteristic volume ~ d\ contributing to the integral in Eq. (12). Thus the 
two effects together vary as (R/d\) 2 d\ oc 1/d 2 , and have no net dependence on R. 

The damping of the cantilever oscillations can be described alternatively as follows. 
It can be seen from Eq. (11) that there exists not only a ^-component, but also an 
x-component of the field E', i.e., a component parallel to the direction of tip motion 
that reduces the cantilever velocity. The loss W of the kinetic energy per unit time 
is proportional to a friction force and the velocity dX(t)/dt: 




(12) 



which, when used with Eq. (11), becomes 



b(qX uj) 2 _ b{VX uj) 2 
Snjldl ~ 2 7 nfid 2 



(13) 



W = -bqE' x (x = 0,z = d u t)dX(t)/dt , 



(14) 
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where the overline denotes an average over a period of the oscillations. The explicit 
form of W can be obtained using Eq. (11), and it is easily verified that W obtained 
in this way is exactly the same as that given by Eq. (13), i.e., W = Wj. 
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FIG. 2. Cantilever-sample geometry for the case of finite sample thickness. 

In the case of finite sample thickness d s we assume that the Debye layer thickness 
Rd is negligible. (For good metals, this thickness is of order 1 A.) For a finite sample 
thickness we must account for the additional boundary condition at z — —d s (see 
Fig. 2). To derive this boundary condition, we use the fact that when d s y Rd the 
static field does not penetrate to the depth d s . An electric charge can be induced at 
the second boundary of the sample only by virtue of the field E' oc uj. Therefore the 
term dp/dt in the continuity equation can be neglected, as it is of order u 2 , so that 
the boundary condition at z — —d s for zero surface charge is 

j z (x, z = -d s ) = (iE z (x, z = -d 8 ) = . (15) 

Thus the solution of the reduced continuity equation V • /xE' = with boundary 
conditions (8) and (15) determines the field E' in the bulk of the sample. Further 
analysis becomes more convenient if we consider the field E' to be induced by charges 
on the surface z = and on the "image plane" at z — —2d s . Setting the charge 
density at the image plane equal to that at the z = plane, i.e., a'(x,z = 0) = 
a'(x,z = —2d s ), ensures that the boundary condition j z (x,z = —d s ) = is satisfied. 
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The electric field E' can be expressed in terms of a' as 

d 



/OO (j 
dx'a'(x',t) — [mlr-e^'l +\n\r-e x x' + 2d s e z \] . (16) 



To obtain an expression for a' we use Eqs. (8) and (16): 

1 , v ■ d ( 1 \ 
qdiuoXs sin oot— 



7r/i^ 1 u <9x ^rr 2 + d\) 

2d s - A A 



2 y dx'a\x',t) 



17) 



_ (24 - A) 2 + (x - x'f A 2 + (x - x'f 
The solution of this equation for a'(x,t), for the practically important case of a thin 
sample (d s <C d±), is 

= -^osin^ 47r2M( ^ + ^ . (18) 

The field component responsible for friction, E' x {x = 0,z = di,t), then follows from 
Eqs. (16) and (18): 

, . . QojXq sin u)t , . 

^^"•^^'^S^r- (19) 

This is seen to be greater by a factor of 2di/d s than the corresponding field for the 
infinitely thick sample. Therefore the Joule loss Wf s is increased by the same factor: 

3 2 ii /2 vr/irf 3 /2rfs • 

The increase of Joule dissipation in a thin substrate is due to the decrease of the 
characteristic volume where recharging occurs. In this case the effective resistivity of 
the substrate is increased, giving rise to additional ohmic dissipation. 



III. SPHERICAL AND ELLIPSOIDAL TIPS 



We now consider the case of a spherical or ellipsoidal tip. As before we assume a 
tip-plane separation much smaller than the curvature radius R, allowing us to make 
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a "proximity approximation" as follows. For two parallel surfaces (the limit R — > oo) 
the electric field has only a z component, and the attractive force between any pair 
of opposing elements is 

A/i = Asif(d) = Asi-^p . (21) 

For large (but finite) R we may assume that Eq. (21) holds, but with d replaced 
by the distance Di between opposing elements Asj (see Figure 3). It follows from 
elementary geometrical considerations that for small angles ctj the distance is given 
by 

(#/2K 2 , (22) 

and the total force acting on the tip is 

™ a , V 2 R 2 r daa , . 

r =£A /|M _/_. (23) 

Since small angles a ~ (d/R) 1 ^ 2 make the dominant contribution to the integral, the 
integration range can be taken from to oo. Thus we have 

RV 2 

T = — . (24) 

It is worth mentioning that the application of this "proximity approximation" ap- 
proach to the case of the cylindrical tip gives exactly the force of attraction (5). 
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FIG. 3. In the "proximity approximation" the force is obtained by adding the pairwise 
forces Afi between opposing surface elements separated by distance D\. 

Let us consider now an ellipsoidal tip surface having two different curvatures in 
perpendicular directions. In the vicinity of a sample the equation describing this 
surface involves a quadratic dependence on x and y: 

Z ^ y) = & + &> (25) 
zn x ZKy 

where R x and R y are the radii of curvature in the two directions. Our previous 
considerations may be repeated with the modification that the tip-surface distance 
Di is now taken to be 

D i = d + z(x i ,y i ) • (26) 

For the purpose of integration it is useful to introduce the change of variables 

x = {R x /R y ) l l 2 x' , y = {R y /R x ) l ' 2 y' , (27) 
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in terms of which the problem is effectively reduced to the spherical case, and the 
force for the ellipsoidal tip is given by Eq. (24) with R = (R x R y ) l l 2 . 

The calculations of Joule loss in the case of a spherical tip proceed as outlined 
in the previous section. The surface charge density of the sample for the moving 
cantilever now depends on both x and y: 

The additional surface charge a' that generates dissipative currents in the sample is 

and the electric field responsible for dissipation is determined from a' by 

E'(r,i) = -v/dr ± ^i, (30) 
where r = (x,y,z) and rj_ = (x,y). The expression for the Joule loss is 

= <w*v» 

J 32 v / 2tt : V 3 / 2 V ; 

where the (dimensionless) integral / is given by 

1 = J J dr ± dv' ± \r ± - r'J" 1 J^(l + r l)' 1 ^( 1 + r 'D~ l « 48 ■ ( 32 ) 

The d -3 / 2 dependence of Wj on d given by Eq. (31) is somewhat weaker than in 
the cylindrical case, and is comparable to the dependence on d of the frictional force 
observed by Stipe, et al. [1]. However, the numerical value given by Eq. (31) is too 
small to explain the experimental results. This question will be taken up in Section 
6. 



13 



IV. CASIMIR FORCE 



The "proximity approximation" employed in the previous section can be applied 
to the case of the Casimir attraction. We consider the case of small separations d 
such that d <C c/u p ,hc/kBT, where u p is the plasma frequency; u p ~ 10 16 s _1 for 
good metals, so that we require d <C 30 nm. For these small separations we can use 
the Lifshitz formula [5] for the force per unit area: 

where e is the dielectric permittivity, assumed to be the same for both plates. The 
explicit form of e in the Drude approximation is 

£ «) = 1 + ^ , (34) 

where v is the electron collision frequency (~ 10 13 s _1 at T = 300K), which is much 
smaller than cu p , and can be neglected when integrating Eq. (33). It follows from 
Eqs. (33) and (34) that 

Using (35), (21), and (23), we obtain 

r <^J^ < 36 > 

for a spherical tip of radius R. In the case of a cylindrical tip we have 

3bRV 2 huj p 

T{d) = 2*dW ' ( } 

This differs from the force obtained for the case of a cylinder in Reference [14] (T(d) oc 
d-(n-i)^ n = 3 or 4). The latter result was obtained using the unjustified assumption 
that the Casimir force satisfies Laplace's equation in the gap between the metals. It 
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is seen from Eqs. (36) and (37) that the dependence of the Casimir force on d is 
stronger than in the electrostatic case. At small values of d the Casimir force can be 
dominant. 



V. ATTRACTION DUE TO SPATIAL FLUCTUATIONS OF 
ELECTROSTATIC CHARGE 

In this section we consider a substrate with a stationary, inhomogeneous distribu- 
tion of the electric charge, which can occur in good dielectrics. Such a situation was 
approximated experimentally [1] by employing a fused silica sample irradiated with 
7 rays. In the course of irradiation, positively charged centers (Si dangling bonds) 
are generated. Randomly distributed positive charges are embedded in a negative 
background, and on average the sample is electrically neutral. We model the sample 
as consisting of microscopically small volume elements AV{. Each element is chosen 
sufficiently small that not more than one charge center is present in it. The electric 
charge qi of each element can be described by 



where = 1 or for an occupied or unoccupied charge center element, respectively, 
and e& is the background charge in the volume AVi. It follows from Eq. (38) that the 
condition of charge neutrality, (g,) = 0, is satisfied when 



where n = (rij), n < 1. In the following we will consider the fluctuations of charges 
in different volume elements i,j to be statistically independent, so that {5qi5gj) = 
for i 7^ j. The mean square of charge fluctuations within a given volume element is 
given by 



qi = n { e + (1 - ni)e b , 



(38) 



e(rn) 



—en , 



(39) 




15 



(5qMi) = (n 2 t )e 2 + 2ee b {n % {\ - m)) + ((1 - n t ) 2 )e 2 . (40) 

It follows from the definition of rij that n 2 = (1 — n^) 2 = 1 — n,, and nj(l — n^) = 0. 
Then using the assumption that n <C 1, we obtain from Eq. (40) the mean square 

(SqiSqt) « ne 2 , (41) 

which shows that a random charge distribution is consistent with Poisson statistics. 
Each individual charge qi has its own image — qi in the bulk of the cantilever tip. Then, 
in the absence of cross terms, an average tip-sample attractive force is determined 
by adding the contributions from all the charge-image pairs. Thus the electrostatic 
attraction A/j between opposing elements AVi is given by 

TIG 2 

Afl = 4[D(x»Vi)-Zi] 2 ' (42) 

where the coordinates x iy y,, z^ give the position of the ith volume element in the 
substrate, and the total force is obtained by summing over all the elements: 

T=^Y l - . (43) 

4 ^ [D(xi, yi) — Zi] 2 1 1 

After replacing the sum by an integral (nXI — > cj d 3 r, where c is the number of 

charge centers per unit volume), we obtain, after integration over z, 

ce 2 r r dxdy irce 2 R 2 r daa 
= ~J J D(^j = ~^2dT 1 1 + ' { ' 

The integral on the right-hand side is logarithmically divergent. Hence we should 
confine the integration to values of a up to a max ~ Ru P /2R, where Ru p is the char- 
acteristic size of the bottom part of the tip. The value of the integral is not very 
sensitive to Ru p , so that we may set R = Ru p - Then the force of attraction is given 
by 



T = -7rce 2 R\n 
2 



1+ 8d. 



(45) 



16 



To illustrate the importance of this force, we put in Eq. (45) the typical values 
c = 10 18 cm~ 3 (as in [1,25], see also references therein) for the defect concentration, 
and R ~ 1 /im, d — 20 nm. In this case the ratio of the force (45) to the Casimir 
force (36) is 10. This ratio decreases with decreasing d. 



x 




T 

FIG. 4. Transverse displacement X(z, t) of the cantilever centerline. Restoring force X'T 
acts on the tip when a / 0. 

Our considerations in this section have ignored the screening of the electric field 
in the dielectric substrate. This can be justified in the case of very small tip-sample 
separations (substantially smaller than screening length) , as only defects in the surface 
layer of thickness d contribute to the integral in Eq. (44). When the screening is 
important, the effective electric field inducing image charge in the metal is decreased 
by the factor of (e + l)/2 (see Reference [24]). This factor is equal to 2.5 in the 
case of silica. In addition, the tip image charge induces another image charge in the 
dielectric, and so forth. Nevertheless, for rough estimates it is sufficient to consider 
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only the first pair of charges, equal to ±2e/(e + 1). 

VI. INFLUENCE OF ATTRACTIVE FORCE ON CANTILEVER 

EIGENFREQUENCIES 

When the cantilever geometry satisfies the conditions L > 6 > a we can employ 
the equation of elasticity [27] to describe the dynamics. In the case of an isotropic 
material, the transverse displacement X(z, t) of the beam centerline (along the x axis, 
see Figure 4) satisfies the differential equation 

Pc S— X(z, t) = -EI t —X(z, t) , (46) 

where p c is the density of the cantilever material, S = ab is the cross-section area, 
E is Young's modulus, and 7j = a 3 b/12 is the bending moment of inertia. The 
clamped end, at z — 0, imposes the boundary conditions X(z = 0) = X'(z = 0) = 0. 
The free end at z = L imposes the boundary conditions M y = ELiX" = and 
F x = —EliX'" = 0, where M y and F x are the moment of the elastic force and the 
elastic force in the y and x directions, respectively. The fundamental eigenfunction 
is given by X{z,t) = X (z) cos(wt), where 

Xq(z) = A[(coskL + cosh kL) (cos kz — cosh kz) + (sin/tL — sinh kL) (sin kz — sinhrez)] . 

(47) 

Here A is a constant, k 4 = u 2 p c S/EIi, and 

ou={Z.™/L 2 ){Eh/p c S) 1 l 2 . (48) 

To estimate u we put a = 250 nm, b = 7 x 10 3 nm, L = 250 pm, and use the 
material constants of Si: p c = 2.33 g cm -3 , E = 9 x 10 11 g cm" 1 s~ 2 . Then we have 



18 



(jj/2tt 4 x 10 3 s 1 , which is approximately the same as the value quoted in Reference 
[!]■ 

In the presence of a stretching force T, Eq. (46) is modified to [27] 

Pc S—X(z, t) = -EI t —X(z, t) + T— X(z, t) . (49) 

Furthermore the boundary condition at z = L (X'" = 0) should be replaced by the 
more general condition —EIiX'" + TX' = 0. The term TX' describes the effect of 
the restoring force when there is a nonvanishing angle between T and the centerline 
direction at z = L (see Figure 4). This value of the restoring force is valid in the case 
of small angles a, when sin a ~ tana = —X'(L). 
The eigenfunctions of Eq. (49) are determined by 

X (z) = A[(A 2 cosA + A 2 cosh A) (cos \z - coshA^) 

+ (A sin A - A sinh A) (A sin \z - A sinh Az)], (50) 

where 

A = L[(r/2) 2 + {k a + r 2 /4)V 2 ]V 2 , A = L[(-r/2) 2 + (k 4 + r 2 /4)V 2 ]V 2 , r = T / EI t . 

(51) 

For simplicity the coordinate z in Eq. (50) is given in units of L. The value of k in 
Eq. (51) obeys the following eigenvalue equation: 

2k, 4 + (r 2 + 2k 4 ) cos A cosh A + k 2 t sin A sinh A = . (52) 

The minimum value of k (K min ) that obeys Eq. (52) determines the fundamental 
eigenfrequency and the explicit form of the fundamental eigenfunction (50). 

In the general case the fundamental frequency is a function of the tip-sample sep- 
aration, bias voltage V, tip geometry, and concentration of charge centers via the 
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attraction force T(V). In the case of attraction due only to an external bias voltage, 
there is a quadratic dependence of T on V. The proportionality coefficient depends 
on the geometry of the tip-sample system. To illustrate the oj(V) dependence we use 
the parameters given at the beginning of this section. The results of numerical calcu- 
lations are shown in Figure 5 for a cylindrical geometry and for the same parameters 
assumed as at the beginning of this section. It is seen that the bias voltage can have 
a considerable effect on the frequency shift. 



co(V) /co(O) 




1< *>.o 

Volts 

FIG. 5. Numerical results for the cantilever frequency to and effective mass m e // as a 
function of the bias voltage V in the case of a cylindrical tip (see text). 

In the limit of very large V (when TL 2 /(EIi) ^> 1), the first term on the right- 
hand side of Eq. (49) can be omitted, and Eq. (49) is transformed into the equation 
for the vibration of a string. In this linear dependence of uo on V (ou oc V/L) 

should occur. 

Along with the resonant frequency, the cantilever effective mass also depends on 
V . The effective mass m e ff is the coefficient entering the oscillator equation 

d 2 x(t) , ax(t) 



+ r- 



+ kX(t) = 



(53) 



as a parameter that determines the inertial force, (r is a friction coefficient, k a 
spring constant, and X(t), as previously, is the coordinate of the tip.) The mass 
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m e f f is related to the coefficients in Eq. (49). This relation can be obtained from the 
requirement that the kinetic energy of the flexible beam equals that of the oscillator, 
the two being one and the same physical quantity. This condition defines m e ff as 



where X (z) is given by Eq. (50). 

Thus m e ff depends on the variation of X with z, which in turn depends on the 
stretching force T and the other parameters in Eq. (49). In the simplest case T = 0, 
we have m e ff = p c LS/4:, i.e., a quarter of the beam's actual mass. In general the 
possible variation of m e ff should be taken into account when the friction coefficient 
T is obtained experimentally from the relation T = m e ff(v/Q, where Q is the quality 
factor. In other words, not only the quality factor and the eigenfrequency, but also 
the value of m e ff is required to obtain T. The dependence m e ff(V) is shown in Figure 
5. It is seen to be similar to that of uo(V) at small values of V. 

Figure 6 shows the effect of the Casimir force on the fundamental frequency in the 
case of a spherical tip. The three curves illustrate the tendency of u{d, R) to increase 
with the attractive force, which varies with d and R as T ~ R/d 2 (see Eq. (36)). The 
frequency shift observed in Reference [1] at d = 2 nm for a gold sample corresponds to 
the curve R = 0.5 //m in Figure 6. When we use the experimental value 1 /im of the 
radius of curvature in the y direction, and set the radius of curvature in the x direction 
equal to the cantilever thickness a = 0.25yum, we obtain R exp = (R x R y ) 1 / 2 = 0.5 /xm., 
i.e., R exp = R. Hence the frequency shift observed by Stipe et at. [1] (see Figure 2c 
of Reference [1]) may be attributable to the Casimir effect. 




(54) 
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Separation d (nm) 

FIG. 6. Effect of Casimir force on the frequency oj. Solid line: R = 0.5 /urn; dashed line: 
R = 1 fim; dotted line: R = 2 fj,m 

It is interesting to compare the experimental and theoretical values of T. The 
Joule loss Wj [see Eqs. (13), (20), and (31)] determines the friction coefficient: It 
follows from energy conservation and Eq. (53) that 

- - 2W > (55) 



In the case of the cylindrical tip we have, from Eqs. (14) and (55), 

bV 2 



(56) 



Putting, as previously, b = 7 x 1CT 4 cm, a conductivity /x = 4 x 10 17 s" 1 for gold at 
300 K, d = 20 nm, and V = 1 Volt, we have T = 2.4 x 10~ 20 kg/s, which is much 
smaller than the experimental value of 3 x 10~ 12 kg/s. For the more realistic case of 
a spherical tip with R = 1 nm, our estimate for T using Eq. (31) is two orders of 
magnitude smaller than that given by (56). 
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VII. CONCLUSIONS 



The analyses presented in this paper indicate that both electrostatic forces and 
Casimir forces can have a strong effect on the cantilever vibrations. The cantilever 
eigenfrequencies depend on the attractive force T, which is very sensitive to the tip 
geometry for small values of the tip-sample separations d. The effect of various forces 
on the cantilever eigenfrequencies appears to us to be an important consideration for 
the practical utilization of cantilever-based devices. It has served as a motivation in 
this paper to study the tip-sample interaction in some detail. 

We have shown that, for small separations, T depends on the shape of the can- 
tilever tip. The dependence of T on the separation d is different for cylindrical and 
spherical tips. Thus, in the case of the electrostatic force due to a bias voltage, the 
attractive force varies as d~ 3 / 2 and d^ 1 for cylindrical and spherical tips, respectively. 
In the case of the Casimir force the dependence of T on d was found to be d~ 5 ^ 2 
and d~ 2 , respectively, for the cylindrical and spherical tips. We have shown that the 
Casimir force is responsible for the frequency shift, which is of the same order as that 
obtained experimentally in [1] 4.5 %) for a gold sample at small separation d — 2 
nm. 

The attractive force depends furthermore on the radii of curvature, and our anal- 
ysis allows for the possibility that the radii may be different in different directions. 
In principle, any variations of tip geometry (for instance, due to adsorption of new 
molecules or blunting after contact with a sample) may be detected by measurements 
of the frequency shift, and our theory, which connects the variations of the attractive 
force with the frequency shift, could be helpful in estimating the character and scale 
of these variations. 

We have calculated the attractive force between a system of randomly distributed 
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positive charges, embedded in a negatively charged background, and a metal tip. In 
the case of an electrically neutral system, only the attraction between each charge 
and its image contributes to the total force. This is a spatially fluctuating interaction 
because the overall force is linear (not quadratic) in the concentration of elementary 
charges. Our estimates, based on the assumption of an uncorrelated distribution of 
charge centers in the bulk of the substrate, indicate that the fluctuating force exceeds 
the Casimir force for parameters appropriate to the experiments reported in Reference 
[1], which employed irradiated silica as a substrate. 

We have derived analytical expressions for Joule losses in a metal substrate. The 
Joule mechanism does not explain the cantilever damping measured in [1], but our 
analysis may nevertheless be useful for other systems with high-resistivity substrates. 
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